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The characterization of an infinite-order quantum phase transition (QPT) by entanglement mea- 
sures is analyzed. To this aim, we consider two closely related solvable spin-1/2 chains, namely, the 
Ashkin- Teller and the staggered XXZ models. These systems display a distinct pattern of eigen- 
states but exhibit the same thermodynamics, i.e. the same energy spectrum. By performing exact 
diagonalization, we investigate the behavior of pairwise and block entanglement in the ground state 
of both models. In contrast with the XXZ chain, we show that pairwise entanglement fails in the 
characterization of the infinite-order QPT in the Ashkin- Teller model, although it can be achieved 
by analyzing the distance of the pair state from the separability boundary. Concerning block entan- 
glement, we show that both XXZ and Ashkin- Teller models exhibit identical von Neumann entropies 
as long as a suitable choice of blocks is performed. Entanglement entropy is then shown to be able to 
identify the quantum phase diagram, even though its local extremes (either maximum or minimum) 
may also appear in the absence of any infinite-order QPT. 



PACS numbers: 03.65.Ud, 03.67.Mn, 75.10.Jm 



I. INTRODUCTION 



The behavior of entanglement in many-body systems 
has attracted great attention in recent years due to its 
promising potential to realizing quantum information 
tasks [l|-[3| as well as its relationship with quantum crit- 
ical phenomena Q. In particular, it has been observed 
that entanglement can identify and characterize a quan- 
tum phase transition (QPT) [1, @. QPTs are associated 
with critical changes in the ground state of a quantum 
system due to level crossings in its energy spectrum, oc- 
curring at low temperatures T (effectively T = 0). Specif- 
ically, a first-order QPT is characterized by a finite dis- 
continuity in the first derivative of the ground state en- 
ergy. Similarly, a second-order QPT - or a continuous 
QPT - is characterized by the existence of an infinite cor- 
relation length and a power-decay of correlations, which 
is often manifested by a finite discontinuity or divergence 
in the second derivative of the ground state energy, as- 
suming the first derivative is continuous. A more sub- 
tle class is the so-called infinite-order QPTs, for which 
any finite-order derivative of the ground state energy is a 
continuous function of the relevant parameters. Promi- 
nent examples of such QPTs are provided by the metal- 
insulator point in the fermionic Hubbard model and the 
SU(2) (antifcrromagnetic) Heisenberg point in the XXZ 
spin-1/2 chain. 

In recent years, it has been noticed that the behavior 
of entanglement (or its derivatives) in the ground state of 
a many-body system undergoing a first-order or a con- 
tinuous QPT exhibits a non-analiticity. Indeed, in the 
case of first-order QPTs, discontinuities in the ground 
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state entanglement were shown to detect the QPT 0-[i- 
For the case of second-order QPTs, the critical point is 
found to be associated with a singularity in the derivative 
of the ground state entanglement, as first illustrated for 
the transverse field Ising chain in R ef. |j| , and generalized 
in Refs. 10-12] (see also Refs. |13l - [l7| for an analysis in 
terms of other entanglement measures). The behavior of 
entanglement for first-order and second-order QPTs have 
also been discussed in general grounds in Refs. [isl. [l9|. 
For infinite-order QPTs, although no independent-model 
analysis is available, entanglement has been found to ex- 
hibit a local extreme (either maximum or minimum) at 
the quantum critical point (QCP). This has indeed been 
shown for the spin-1/2 XXZ chain [20. - 2^^] and for the 
fermionic Hubbard model [131 ■ Whether or not this is a 
general property of an entanglement measure for a con- 
venient partition of the system remains unresolved. 

In order to make progress on this matter and on 
the general properties of an infinite-order QPT, we will 
consider in this paper the characterization of entangle- 
ment into two closely related solvable spin-1/2 chains, 
namely, the Ashkin- Teller and the staggered XXZ mod- 
els. The Ashkin- Teller model has been introduced as 
a generalization of the Ising spin-1/2 model to investi- 
gate the statitiscs of two-dimensional lattices with four- 
state interacting sites [Slj . Since then, it has attracted a 
great deal of attention due to its motivations in a wide 
range of research fields. First, both its classical and 
quantum versions exhibit a rich phase diagram [32l . [33| , 
which makes the model a prototype for the investiga- 
tion of phase transitions and critical phenomena. For 
instance, it has recently been shown that the quantum 
Ashkin- Teller model in one-dimension exhibits an exam- 
ple of disorder rounding of a first order quantum phase 
transition into a continuous phase transition [34] . which 
may find applications to numerous complex strongly cor- 
related systems. Second, it has been experimentally 
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realized by magnetic compounds formed by layers of 
atoms adsorbed on clean surfaces, e.g., Selenium ad- 
sorbed on Ni(lOO) surface (sH]. Moreover, the univer- 
sality properties Aslikin- Teller model may also be re- 
lated to many other intereseting applications, such as 
nonabelian anyons models (36l |. orbital current loops in 
Cu02-plaquettes of high-Tc cuprates [s^l , and elastic re- 
sponse of DNA molecules [S^ . Remarkably, the Ashkin- 
Teller model and the XXZ chain display a distinct pat- 
tern of eigenstates but exhibit the same thermodynamics, 
i.e. the same energy spectrum. In particular, the same 
quantum phase diagram applies to both models, although 
entanglement can be distinct in each case. Our aim is 
then analyze this quantum phase diagram and investi- 
gate to what extent entanglement measures will agree on 
the description of the infinite-order QCP. Indeed, as we 
will see, in contrast with the XXZ chain, pairwise entan- 
glement between nearest neighbors in the Ashkin- Teller 
model is not able to identify the QPT. By focusing on 
block entanglement, we will then show how to reconcile 
the behavior of entanglement in both cases. Moreover, 
we will show that local extremes may also occur in the 
absence of any infinite-order QPT, indicating that fur- 
ther analysis of the critical behavior is demanded in a 
scenario where entanglement displays a local maximum 
or minimum driven by a relevant parameter. 



II. THE ASHKIN-TELLER MODEL AND ITS 
MAP INTO THE STAGGERED XXZ CHAIN 

Let us begin by introducing the quantum Ashkin- Teller 
model in one-dimension, whose Hamiltonian for a chain 
with M sites is given by 

M 

M 

where cr" and r" {a = x,y,z) are independent Pauli 
spin-1/2 operators, J is the exchange coupling constant, 
A and /3 are (dimcnsionless) parameters, and periodic 
boundary conditions (PBC) are adopted, i.e., cr^z+i = 
and r^_|_i = {a — x,y,z). The Ashkin- Teller model 
is Z2 Z2 symmetric, with the Hamiltonian commuting 
with the parity operators 

M M 

Vi = l[<7^ and ■P2 = l[r^- (2) 

Therefore, the eigenspace of Hat can be decomposed 
into four disjoint sectors labelled by the eigenvalues of 
Pi and P2, namely, Q = {Vi = -Hl,p2 = +1), Q = 1 
(Pi = +l^V2 - -1), Q = 2iVi^ -l,V2 = -1), and 
Q = 3 CPi = -l,'P2 = +1). By the symmetry of Hat 



under the interchange cr" O r", the sectors Q = \ and 
(5 = 3 are degenerate. Moreover, we observe that the 
ground state belongs to the sector (5 = 0. 

In order to map the Ashkin- Teller model into the stag- 
gered XXZ chain, we consider two sets of 2Af link vari- 
ables 7j|j = 1, ■ ■ ■ , 2Af}, which are defined by 

?72j-l =crj, 72j-i=t/, 

?72j = Cr|cr|+i, 72j = T^T^+l- (3) 

These variables satisfy the conditions 

7y| = l, 72 = 1 (j = l,...,2Af). (4) 

Moreover, they obey the algebra 

[Vj^lk] = 0, [77j,?7fc] = 0, [7j,7fc] = 0, (5) 

for |j - fc| ^ 1 and (j, k) ^ (1, 2M) and (j, k) ^ (2Af, 1), 
while 

['7i,7fc]=0, {??i,%} = 0, {7j,7fc} = 0, (6) 

for |j ^k\ = l or [j,k) = (1,2A/) or (j, fc) = (2Af,l). 
However, note that, since PBC are adopted, the link 
variables are not completely independent, turning out to 
obey the constraints 

M M 

X{mi^X{i2, = t. (7) 

Moreover, by writing out Hat in terms of the link vari- 
ables, we obtain 

2M 

Hat ^ ^'^J (^2j-i + 72i-i + A772j_i72j-i) 

2M 

-JPY. ('?2j + 12, + Ar72j72j) ■ (8) 

The equivalence between the Ashkin- Teller model and 
the XXZ chain can be established by considering the stag- 
gered spin-1/2 XXZ chain with 2Af sites, whose Hamil- 
tonian reads 

M 

HXXZ = -Y.'^ K-l4j + ^1,-1^1, - A^72,_1^2j] 
M 

-J/3 + 44+1 - , (9) 

where, as before, J is the exchange coupling constant, 
A and /? are the anisotropy and staggering (dimension- 
less) parameters, respectively, and PBC are adopted, i.e., 
<J2M+i = (a = X, y, z). Note that, differently from the 
Ashkin- Teller model, the XXZ chain has one spin-1/2 
particle per site. The XXZ model is U(l) invariant, with 
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the Hamiltonian commuting with the total spin operator 

[Hxxz,S']=Q. (10) 

Therefore, the eigenspace of Hxxz can be decomposed 
into 2M + 1 disjoint sectors labelled by their correspond- 
ing magnetization quantum number n — M — r, with 
r ~ 0,1, ...,2M denoting the number of spins reversed 
from the state with all spins down. We observe that the 
ground state belongs to the sector n = 0. A schematic 
view of the XXZ and Ashkin- Teller chains is shown in 

Fig.m 
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FIG. 1: (Color online) (a) Ashkin-Teller spin-1/2 chain. The 
lattice is composed by two independent spin-1/2 particles per 
site j described by Pauli operators {(T",r"}. (b) XXZ spin- 
1/2 chain. The lattice is composed by one spin-1/2 particle 
per site j described by Pauli operators {cr^}. 

As defined for the Ashkin-Teller model, we introduce 
two sets of link variables, {?7j, 7j|j = 1, • • • , 2M}, which 
are given by 

V2j-i = <^2j-i^2j, 72J-1 = 

V2j = (tI- gI-^^ , = cr^j- cr^j-^ 1 , (11) 

with j — 1, • • • , M. These variables also satisfy the con- 
ditions given by Eq. (|4]) and the algebra given by Eqs. ([5]) 
and Moreover, we observe that the constraints 

obeyed now by the link variables with PBC adopted read 

M M 

n '72j-l72i = n T2i-l7?2i = 1. (12) 

By writing Hxxz m terms of {ryj, 7^} we obtain the same 
expression as in Eq. (jS]). Therefore, an equivalence be- 
tween Hxxz and Hat can be achieved as long as the con- 
straints given be Eqs. (O and (IT^ can be made compat- 
ible between each other. In this direction, we can show 
that, despite displaying rather different ground state vec- 
tors, the XXZ and Ashkin-Teller models, as defined by 
Eqs. dU and H]), exhibit the same ground state energy 
More specifically, with PBC adopted, it can be shown [25'| 
that the energies in the sector Q = of Hat: which 
contains the ground state, will occur in the spectrum of 
Hxxz- Indeed, observe first that eigenstates of Hat into 
the sector Q — Q are characterized by 

M M 



which is a consequence of Eq. ([2]). Then, by inserting 
Eq. (fT3| into Eq. ([7]), we exactly obtain the constraint 
given by Eq. (|12p . Therefore, all the energy levels belong- 
ing to Q = also appear in Hxxz- Conversely, we can 
show here that the ground state energy of the XXZ model 
is also contained into the spectrum of Hat- Indeed, be- 
sides the U(l) symmetry given by Eq. (fTOl) . Hxxz is also 
invariant under Z{2) transformations, namely, it com- 
mutes with the parity operators 

2M 

with a = x,y, z. In particular, the ground state of Hxxz 
is into the sector of eigenstates such that 

Q,|^)=-f|^), Q,|V) = +|V'). (15) 

Then, by using Eq. (IT5|l into Eq. p^ . we exactly obtain 
the constraint given by Eq. ([7]). Therefore, all the en- 
ergy levels belonging to the sector Q.j^ — -1-1, = -1-1 
also appear in Hat- Hence, we can conclude that the 
ground state energy Eo{A,P) of both models are iden- 
tical for any A and /3. Remarkably, this equivalence 
can be extended to the whole spectrum if open ends are 
adopted [2^. Bearing in mind that Hxxz and Hat 
have the same ground state energy, they will manifest 
the same quantum phase diagram as the anisotropy A 
and the staggering parameter /3 are varied. Naturally, 
the character of the quantum phase itself depends on the 
particular model, since it is associated with the prop- 
erties of the ground state vector. The quantum phase 
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FIG. 2: (Color online) Quantum phase diagram of Ashkin- 
Teller and staggered XXZ models. The energy scale J is set 
to one. 

diagram for Hat and Hxxz is sketched in Fig. [2] The 
quantum phases for the Ashkin-Teller model can be de- 
scribed as the following: the ferromagnetic phase (F), 
where = (rj) 7^ 0, the paramagnetic phase (P), 

where ((t|) = (rj) = 0, and the partially ordered phase 
(G), where {<^jTj) 7^ 0. Since our main focus is the in- 
vestigation of the infinite-order point /3 = 1 and A = 1, 
we will concentrate on the diagram starting on point 
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A = — -\/2/2, where we have the beggining of a critical 
hne of continuous QPTs. 



III. PAIRWISE ENTANGLEMENT 

Let us first analyze the behavior of entanglement for 
pairs of spins in the XXZ and Ashkin- Teller models. To 
this aim, we quantify entanglement by employing the 
negativity [1^ |23| , which is given by 



U{p'^) = 2 max(0, - min(A^^)), 



(16) 



where are the eigenvalues of the partial transpose 
pi-j^TA Qf ^Yie two-spin density operator p*-', defined as 
(a/3 1 p^* 1 7(5) = (7/3 1 p\a5). For the (non-staggered) XXZ 
model, pairwise entanglement between nearest-neighbors 
in the ground state of the chain has been evaluated in 
several previous works [20l - |23j . with entan glem ent mea- 
sured either by negativity or concurrence [23] ■ For the 
specific case of ground state entanglement in the XXZ 
model, negativity and concurrence turn out to be iden- 
tical [23, [231 ■ In particular, at the infinite-order QCP 
A = 1, it has been shown that the negativity achieves 
a maximum given by 0.386. Indeed, this maximum is 
rather robust, appearing not only in the ground state 
but also for all conformal towers associated with Hxxz 
formed by an infinite number of excited states [23{ . 

Concerning the Ashkin- Teller model, we have two dif- 
ferent spins at each site, given by the Pauli operators a 
and r. We will consider here entanglement between pairs 
Gj — Tj+i, Oj — CTj+i, and r, — Tj+i at nearest neighbour 
sites as well as pairs aj — Tj at the same site (frontal 
pairs). Since the Hamiltonian is symmetric by an inter- 
change cr r, both Uj — (Jj+i and Tj — Tj+i exhibit the 
same entanglement properties. The negativity for these 
pairs is rather different of that for nearest neighbor spins 
in the XXZ chain, displaying no signature (e.g., a max- 
imum) of the infinite-order QPT. This result is plotted 
in Fig. ISja). For the spins aj — Tj and aj — Tj+i, we 
can show that the negativity is vanishing for any A, as 
given by Figs.jH^b) andjU^c). Therefore, entanglement for 
these pairs do not exhibit any indication of the quantum 
critical behavior. However, an indication of the infinite- 
order QPT can be obtained by looking at a closely related 
quantity, which is the distance of the separability bound- 
ary (DSB), denoted by Aip'^). By following Ref. [30|], we 
can define A{p^^) from Eq. through 



-2min(A*i; 



(17) 



Note that A > implies an entangled state with A(p*-' ) = 
Af{p^^). It can be shown that A = for pure separable 
states while A < for mixed separable states. Results for 
A in a chain with 20 spins are shown in Fig. [Sja), [SKb), 
and 12^0) . Note that the distance A for the frontal pair 
of spins presents a cusp at the infinite-order QCP. This 
cusp is actually hidden by the max operation and is kept 





FIG. 3: (Color online) Negativity Af (solid black line) and 
DSB A for the pairs (a) aj — (Jj+i and r, — r^+i , (b) aj — Tj+i , 
and (c) aj — Tj in the Ashkin- Teller chain with 20 spins. Note 
that the only nearest spin pairs with nonvanishing entangle- 
ment are given by curve (a). The results are plotted for /? = 1. 
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FIG. 4: (Color online) Negativity JV (solid black line) and 
DSB A for a frontal pair aj — Tj for chains with 6, 8, and 20 
spins (from top to bottom) as a function of A. Note that the 
operation max hides the nonanalyticity exhibited by A. The 
results are plotted for /J = 1. 



unchanged for larger chains, as shown in Fig. 21 In order 
to understand the nonanalyticity of A in A = 1, let us 
consider the density matrix for aj — Tj , which is given by 



/ u(A) 
v{A) 

v{A) 





(18) 



im(A) + iG(A),t;(A) = i-iG(A), 
iG(A), with m(A) denot- 



1 

" 4 ' 2 
1 - im(A) 



where m(A) 
and w{A) ~ ^ ^ 
ing the magnetization density in the x direction, namely, 
m(A) = {a-} = (rjf) and G (A) = (a^r^). Since the 
density matrix is diagonal, the partial transposition keeps 
the same eigenvalues as p^'^ . Then, it can be numerically 
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shown that the DSB reads 
A = 



i+m(A)-iG(A) 
| + iG(A) 



if A < 1, 
if A > 1. 



(19) 



Therefore, there is a change in the behavior of A exactly 
at A = 1, which makes the DSB a useful pairwise quan- 
tity capable of identifying the infinite-order QCP in the 
Ashkin- Teller model. 



IV. BLOCK ENTANGLEMENT 

As we have seen, in contrast with the XXZ model, 
pairwise entanglement alone (with no DSB suplementary 
analysis) fails to identify the infinite-order QCP in the 
Ashkin- Teller model. In this section, we will show that 
bipartite block entanglement as measured by the von 
Neumann entropy is able to provide a unique descrip- 
tion of entanglement for both the XXZ and Ashkin- Teller 
models. Given a quantum system in a pure state \'ip) and 
a bipartition of the system into two blocks A and B, en- 
tanglement between A and B can be measured by the 
von Neumann entropy S of the reduced density matrix 
of either of blocks, i.e., 



S = -Tr {pA log2 Pa) = -Tr (ps log2 ps) 



(20) 



where pA — Tr^/O and ps = Tr^p denote the reduced 
density matrices of blocks A and B, respectively, with 

Concerning the XXZ model, it has been shown that 
local extremes of the von Neumann entropy are able to 
identify the quantum phase diagram (2]| . In the case 
of the Ashkin- Teller model, we will show here that the 
entanglement entropy can also identify the infinite-order 
QPT by a local extreme, displaying a superior behavior 
in comparison with pairwise entanglement. Moreover, for 
a suitable choice of blocks, we can also show that the en- 
tanglement entropy is the same as that of the XXZ model. 
Indeed, a sublattice composed by a contiguous block of 
spins in the XXZ chain will exhibit the same von Neu- 
mann entropy as a set of contiguous frontal pairs of spins 
in the Ashkin- Teller model. This is analytically worked 
out in the Appendix for a two-spin block and numeri- 
cally checked for larger blocks. The map of equivalent 
entropies is skechted in Fig. [S] 



• • 

Ashkin-Teller 

FIG. 5: (Color online) Map of spins to obtain the equivalence 
of the entanglement entropies for the Ashkin-Teller and XXZ 
models. 

Note that the map above refers to the von Neumann 
entropy only. Naturally, it is not applicable to obtain 
equivalences in the case of other physical quantities, e.g.. 



pairwise entanglement. As we have seen, no equivalence 
can be obtained in that case. In order to investigate 
the characterization of the QPT by using bipartite block 
entanglement, let us first consider the situation (3 — 1 
(non-staggered model). Then, by performing exact diag- 
onalization, we plot in Fig. [Blthe von Neumann entropy 
S'(A) for a single frontal pair of spins in a chain with 
different lengths as a function of the anisotropy A. Note 
that, independently of the chain size, S'(A) presents a 
maximum exactly at A = 1 (see Inset). We can also show 
that the maximum at the critical point is robust against 
the increase of the size of the blocks as well as the choice 
of the sublattices. This is particularly shown in Fig. [3 
where it is exhibited the entropy in the Ashkin-Teller 
model for different choices of sublattices in a chain with 
20 spins. Note that the entropy gets larger as we increase 
the number of bonds between the sublattices, i.e., ^(A) 
increases according to the direction red [curve (a)] to blue 
[curve (c)] of Fig. [T] Again, a maximum is observed at 
A = 1 (see inset). As mentioned before, entanglement in 
the Ashkin-Teller model will match entanglement in the 
XXZ chain for contiguous blocks of spins (e.g., the red 
plot [curve (a)] of Fig. [7]). 




FIG. 6: (Color online) Block entanglement S'(A) as a function 
of the anisotropy A for a single frontal pair of spins. From 
bottom to top, the chain lenghts are 6, 8, 10, 12, 14, and 20 
spins, respectively. Inset: the derivatives of ^(A) vanish at 
A = 1. 

Hence, the von Neumann entropy provides a com- 
pletely equivalent description of block entanglement in 
the XXZ and Ashkin-Teller models, being able to char- 
acterize the QPT by a maximum at A = 1. We shall 
now analyze the case of staggering parameters (3^1. 
Let us consider the entanglement between a frontal pair 
and all the rest of the chain in the Ashkin-Teller model. 
The resulting entropy for a chain with 20 spins is plot- 
ted in Fig. [8] as a function of A for several values of the 
staggering parameter /?. 

Note that for /3 > 1 the von Neumann entropy dis- 
plays a maximum at A = 1 while for /3 < 1 it is char- 
acterized by a minimum at A = 1. The change in the 
concavity of S'(A) at /3 = 1 is rather robust, being inde- 
pendent of the size of the chain. However, by looking at 
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ings at the phase diagram for a fixed p. Indeed, these P-G 
and F-G crossings are second-order QPTs, with the first 
derivative of entanglement getting a pronounced maxi- 
mum (or minimum) at the QPT. This is exhibited in the 
Inset of Fig. [8] The behavior above is also kept if we 
increase the size of the block. Indeed, we plot in Fig. [S] 
the entanglement between a quartet composed by two 
frontal spins and all the rest of the chain, which exhibits 
the same pattern of maxima and minima for the entropy 
as in the case of single frontal spin pairs. 



FIG. 7: (Color online) Block entanglement S(A) as a function 
of the anisotropy A for L = 4 sites in different configurations 
in a chain with TV = 20 spins. Inset: derivatives of S'(A), all 
of them vanishing at A = 1. 




FIG. 8: (Color online) Entanglement entropy for a frontal pair 
(jj — Tj in the Ashkin- Teller chain with 20 spins as a function 
of A for different values of the staggering couplings p. From 
bottom to top, the figures are plotted for /3 = i, |, 1, |, 
and |, respectively. Inset: the derivatives of S'(A) vanish at 
A = 1 and indicate a pronounced maximum (or minimum) 
around the second-order QPT. 



the quantum phase diagram at Fig. [51 we can see that, 
for /3 7^ 1, there is no infinite-order QPT at A = 1. 
Therefore, we are providing here an example of a local 
extreme of entanglement that is not associated with a 
QPT. Nontrivial staggering is then able to introduce a 
phenomenon that is absent in the standard case. Never- 
theless, observe that the concavity of the curve can char- 
acterize the paramagnetic- ferromagnetic QPT, namely, 
the von Neumann entropy always displays a maximum 
at the ferromagnetic (F) case while a minimum is found 
in the paramagnetic (P) region. Therefore, the concav- 
ity of 5'(A) provides a necessary condition (that is also 
sufficient for A < 1) to determine the phases F and P. 
Moreover, note that entanglement also detects both the 
P-G and F-G QPTs of Fig. [21 given by horizontal cross- 




FIG. 9: (Color online) Entanglement entropy for a quartet in 
the Ashkin- Teller chain with 20 spins as a function of A for 
different values of the staggering couplings /3. From bottom 
to top, the figures are plotted for p — i, |, 1, |, and |, 
respectively. Inset: the derivatives of S'(A) vanish at A = 1 
and indicate a pronounced maximum (or minimum) around 
the second-order QPT. These maxima are expected to scale 
as we increase the size of the chain. 



We observe that the pronunciation of the maximum 
or the minimum of the derivative of entanglement at a 
second-order QPT is expected to evolve to a nonanalyt- 
icity as we increase the size the chain [l^ [3 ■ This is 
indicated in Fig. [TUl where the entropy as a function of 
P reveals the second-order QPTs by vertical crossings at 
the phase diagram given by Fig. [5] for a fixed A. 

Note from Fig.[TU]that, for any size of the chain, entan- 
glement entropy for a quartet tends to 2 as /? — 7> oo. This 
can be understood as a consequence of the map shown in 
Fig. [51 where we replace the entropy of the quartet by the 
entropy of four contiguous spins in the XXZ chain. In- 
deed, in the limit /3 — > oo, the ground state of the Hamil- 
tonian Q is given by a set decoupled dimers, whose den- 
sity operator is p = p2,i®Pi,b®- ■ ■ P2M,i, with /02j,2j+i = 
(1/4) (l2i «)l2j+i <^o-f^+i <8)a|j+i - trf, 0ct|j+i ). 
With no loss of generalility, by analyzing the quartet 
Pi, 2, 3,4, we obtain that pi, 2,3,4 = (l/4)(li (g) p2,3 ® 14)- 
Then it follows that 5'(pi_2,3,4) = 2, which comes from 
the contribution of the spin singlets (2M, 1) and (4, 5). 
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FIG. 10: (Color online) Entanglement entropy for a quartet 
in the Ashkin- Teller chain with 8 (dashed red curve) and 20 
spins (solid black curve) as a function of /3 for A = 5. Inset: 
the first derivative of S{P) displays pronounced maxima at 
P « 0.337 and /3 « 2.14, which are finite-size precursors of 
the second-order QPTs. These maxima are expected to scale 
as we increase the size of the chain. 

V. CONCLUSIONS 

In summary, we have presented a detailed analysis of 
entanglement into two closely related models, namely, 
the staggered XXZ and Ashkin- Teller models. In partic- 
ular, the behavior of the entanglement properties in both 
models has been considered at the infinite-order QCP, 
where we have shown that pairwise entanglement in the 
Ashkin- Teller model fails in the characterization of the 
QPT (with no DSB supplementary analysis). However, 
the von Neumann entropy can offer a unified description 
of quantum criticality, with this equivalence achieved for 
contiguous blocks of spins. Moreover, we have also shown 
that the pattern of local extremes at the infinite-order 
QCP has been obtained by using bipartite block entropy 
for different choices of sublattices, even though local ex- 
tremes also appear in regions where no QPTs occur. This 
result establishes that the maxima and minima of entan- 
glement entropy per se constitute, in the absence of sup- 
plentary analysis, insufficient conditions to ensure the 
existence of an infinite-order QPT. 

The local extremes parttern of entanglement for 
infinite-order QPTs has been observed for a variety of 
different systems, despite no analytical derivation of this 
behavior is available. This is in contrast with first-order 
and continuous QPTs, where the relationship between 
QPTs and nonanalytical behavior of entanglement has 
been generically understood for finite dimensional bipar- 
tite fis'l and multipartite [s^ systems as well as for con- 



tinuous variable models (4C| . We have provided here new 
indications that local extremes provide a typical but not 
sufficient characterization of infinite-order QPTs, which 
may provide a hint for a future derivation of this phe- 
nomenon in general grounds. 
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Appendix 

Let us show in this appendix that the reduced den- 
sity operator p^'^ for a frontal spin pair at site j in the 
Ashkin- Teller is equivalent to a nearest-neighbor two-spin 
reduced density operator pf^j^f in the XXZ chain. Start- 
ing from the Ashkin- Teller model, p^^ reads 

pf ^ \{t4 + ua^ +UTf+v ajr/) , (21) 

where I4 is the 4-dimensional identity operator, u = 
(ctJ) = (Tj), and v = {a^r^). By rewriting Eq. ([2T|) 
in terms of the link variables given by Eq. (jS]) we obtain 

p/^ = I (I4 + uri2j-i + UJ23-1 + V 1123-1123-1) , (22) 

with u = {ri23-i) = {123-1) and v = {1123-1123-1) ■ Now, 
let us turn to the XXZ chain, whose two-spin reduced 
density operator reads 

Pff+i = I {U+pa^a^^, +V^V,+i + . 

(23) 

where p = (crjcrj+i) = (^jf^j+i)' and q = {u]g]^^ . In 
terms of the link variables given by Eq. (1111) , we obtain 

1 / >. 

Pj = ^ [U+pmj-i +PI23-1 - qv2j-ii2j-i) , 

(24) 

with p = imj-i) = {l2j-i) and q = - {r]2j -1123-1)- 
Hence, since u — p and v = —q, we have that Eq. p2|) is 
equivalent to Eq. (|24|). Therefore, the eigenvalues of 
are identical to the eigenvalues of pf^j^f, which means 
that the von Neumann entropy are the same in both 
cases. 
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